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The content of phase information of an arbitrary phase- 
sensitive measurement is evaluated using the maximum like- 
lihood estimation. The phase distribution is characterized by 
the relative entropy-a nonlinear functional of input quantum 
state. As an explicit example the multiple measurement of 
quadrature operator is interpreted as quantum phase detec- 
tion achieving the ultimate resolution predicted by the Fisher 
information. 



There are many approaches addressing the problem 
of quantum phase measurement nowadays. Besides the 
purely theoretical phase concepts anticipating the ex- 
istence of quantum phase as an observable conjugated 
canonically to the number (or difference number) opera- 
tor, there are several operational treatments addressing 
the problem of phase shift measurement within the quan- 
tum mechanics. Particularly, two methods how to de- 
rive the phase information from the phase sensitive mea- 
surement of quadrature operator have been proposed re- 
cently. The former one, so called "phase (measurement) 
without phase (states)" was formulated by Vogel and 
Schleich Q. The method is motivated by the geomet- 
rical meaning of the quadrature- and ideal phase mea- 
surements in phase space. The quadrature eigenstates 
rotated by an angle are used to define a phase distribu- 
tion of a single mode of the radiation field: A balanced 
homodyne-detection scheme measures the electric field- 
strength (variable x) probability 



p(M) = K#c> e | 



(1) 



in dependence on the actual phase of local oscillator 9. 
This quantum detection may be interpreted as measure- 
ment of rotated quadrature operator 



X{9) = -L[ae- 



(2) 



The probability of finding zero electric field plotted ver- 
sus local oscillator phase 9 



P(6) = \(ij\x = 0) e \ 



(3) 



constitutes the proposed phase distribution on the inter- 
val [0, 7r). The phase sensitive data (|l|) resulting from the 
homodyne detection have been interpreted in different 



way by Beck, Smithey and Raymer [@J. Using the opti- 
cal homodyne tomography method |1, the density ma- 
trix may be reconstructed and represented in the phase 
space. Particularly, the authors used the representation 
by Wigner function W(x,p) and linked the phase distri- 
bution to the marginal distribution of Wigner function 



P(<P) 



r drW(x = r cos ip,p = r sin c 



(4) 



The resulting phase distribution is then periodic on the 
interval [0,27r). Nevertheless, such an approach is suffer- 
ing by formal flaw. Since the "probability distribution" 
(Q) yields negative values for superposition of coherent 
states (so called "Schrodinger cat-like states") j|, the 
corresponding operator measure is not positively defined. 
The procedure cannot be therefore interpreted as any 
generalized measurement ||. To get physically reason- 
able interpretation, another distribution function as for 
example the Q-function should be used. The purpose 
of this Rapid Communication is to evaluate the phase 
information included in the phase-sensitive data using 
maximum likelihood estimation. The phase distribution 
then yields the ultimate resolution corresponding to the 
Fisher information. As an explicit example, the quadra- 
ture measurement is interpreted as quantum phase mea- 
surement. The proposed method deals with the observed 
data (|j) in the optimum way. 

Let us formulate the problem for an arbitrary multiple 
measurement of discrete phase-sensitive observable [pi. 
The case of quantum observable with continuous spec- 
trum will be obtained by a straightforward limiting pro- 
cedure. Assume the quantum measurement of quantum 
variable Y yielding discrete spectrum \y k ) enumerated 
for brevity by a (multi)index k. The purpose of phase 
detection is to determine the non-random c-number dis- 
placement parameter 9 in the given interval entering the 
phase displacement transformation || of quantum state 
as \4>(9)) = e^ ieN \tfj), N being a Hermitian operator. 
The variable 9 represents the true value of the phase shift. 
The estimation on the interval 9 £ [0, 2ir) will be con- 
sidered for concreteness. The probability of finding the 
complex amplitude y k by performing the measurement in 
transformed quantum state \ip{9)) is given by quantum 
mechanics as 

p k (e) = \(i>\e ie *\y k )\ 2 . 

Knowing all these probabilities in dependence on the in- 
duced phase shift, an unknown phase shift may be in- 
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ferred on the basis of multiple output data m, y%, . . . ,y n . 
Following the approach used in the Ref. ||, the condi- 
tional phase distribution of inferring phase shift <fi when 
is true, is given by the normalized likelihood function 
as 



1 



C n (0) 



EM*)] 
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The normalization is C n (9) = d<f> {Y[ k \pk{<t>)] Pk{e) } n 
and index k exhausts all the possible values appearing 
with nonzero probability. The number of samples n is 
assumed to be sufficiently large in order to get statisti- 
cally significant sampling. The likelihood function may 
be expressed using the relative entropy 



P{<j>\6) cx e ~ nS(m . 



(6) 



(7) 



The case of phase sensitive observables with continuous 
spectrum y may be easily incorporated in this step defin- 
ing the relative entropy as 



S(<f>\e) = - / dyp(y,9) \np{y, 



(8) 



The preferred phase shift is given by the true value 9, 
since the relative entropy has minimum at S(<f> = 
due to the Gibbs inequality S 



(9) 



The estimation is biased, but may be sometimes well ap- 
proximated by the Gaussian distribution with the vari- 
ance predicted by the Fisher information. Using the Tay- 
lor decomposition of In pk(4>) at the point cf> = 6 the rel- 
ative entropy (||) reads 

sm « - Y,{pk(0) i*PkV) ~^§P[</> - o] 2 + ...}. 



2 Pk{0) 



(10) 



The prime denotes the derivative p' k (9) — dpk{4>) /d4>\^=e- 
The first term represents the Shannon entropy 

S(0) = -£> fc (0) \n Pk (9), 



whereas the second one is the Fisher information 

[p' k m 2 



'(«) = £■ 



Pk(9) 



The variance of phase distribution in this approximation 
is simply 



= 
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Provided that the Gaussian approximation cannot be 
used, the phase resolution may be always evaluated using 
dispersion 



D(6) = y/l-\(e**)\*. 

Here the averaging over the phase is performed in the 
specified phase interval. Dispersion depends on the true 
value of phase shift, in general. For sharp measurements 
it corresponds to the ordinary notion of variance D ss A</> 
restricted to the finite interval . 

As an explicit example assume now the quantum mea- 
surement of phase-sensitive quadrature component (j^) 
performed for concreteness in the coherent state with 
the complex amplitude a = \a\e i(p ; \ip) = D(a)\0). Here 
the displacement operator is D(a) — exp(aa^ — a*a). 
The phase shift of single-mode field is generated by the 
photon-number operator N = a)d. The probability of 
finding the value x of rotated quadrature operator (||) 
may be specified for the given signal state as 



p(x,6') = -^exp{- 

\/7T 



f?\a\ cos 9'} 2 }, 



(11) 



where 9' = 9 — ip is the phase difference between lo- 
cal oscillator and signal fields. The quantum estimation 
problem is the following: The distribution ( fTl] ) is explic- 
itly known as a function of quadrature phase difference 
9' and quadrature component x. These dependencies are 
always experimentally measurable and for example may 
be scanned in advance. The particular choice of Gaussian 
distribution represents an easy example consistent with 
the assumptions of Refs. |Q and Using this knowl- 
edge, an a priori unknown fixed phase difference should 
be inferred as accurate as possible on the basis of limited 
number of measured data Xi,x%, ■ ■ ■ , x n . The straightfor- 
ward application of the theory yields the relative entropy 
as 

S((j)\9') = -\nir + / dxp{x,6')\x- V2|a|cos0] 2 



(12) 



= iln7r+ i + 2|a| 2 [cos0 - cos6»'] 2 , 
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where (f> is estimated (inferred) phase difference. The 
phase distribution inferred after n trials then reads 

P((j>\9') = C(a, 9') exp{-2n|a| 2 [cos0- cos 6>'] 2 }, (13) 

where the normalization is 

C(a,9') = { / d0cxp{-2n|a| 2 [cos(/.-cos6i'] 2 }}- 1 . 
Jo 

Significantly, the resulting phase distribution is not shift 
invariant depending only on the difference 4> — 9' . In gen- 
eral, the phase estimation is biased. It may be interpreted 
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using the standard definition of the normal distribution 
on the circle (von Mises) || 



/(*) = 



1 



2ttI (k) 



exp[« cos (2 — /?)], 



(14) 



x, /3 £ [0, 2tt), which is centered at (3 and characterized by 
the dispersion D 2 = 1 — [7i(k)//o(k)] 2 ; Iq(k),Ii(k) being 
the modified Bessel functions. The phase distribution 
( |l3| ) may be written as the normalized product of two 
von Mises distributions 



p{4>\e') cc /i(0/ 2 (20 



(15) 



centered at /3i = 0, /?2 = — 7r. The parameters are k± — 
An\a\ 2 cos 8' , K 2 — n\a\ 2 . 

Let us detail the phase information included in such a 
phase measurement. The phase distribution ( [l3"| ) exhibits 
the mirror symmetry since P(<f>\9') — P(2ir — <p\6'). Hence 
the phase measurement yields the one-peak distribution 
on the interval [0, 2tt) only if 9' = or 6' = tt. These two 
possibilities are of course distinguishable by the sign of 
the measured quadrature components the prob- 

ability distribution ( [Tl|) indicates. Unfortunately, the 
phase measurement near the points 9' = or 9' = tt 
yields rather bad resolution, as will be seen in the fol- 
lowing. In all the other cases of phase differences 9', the 
detected phase sensitive data corresponding to the statis- 
tics ([ll]) do not distinguish between the values 9' and 
2n — 9' and therefore neither the inferred phase distri- 
bution ( |l3| ) does. This ambiguity of phase measurement 
may be avoided estimating the phase difference on the 
half-width interval [0, tt) only. Normalization should be 
changed to the half value in comparison to the multi- 
plicative factor in distribution (|l3|). The inferred phase 
distribution in dependence on the true phase shift 9 is 
plotted in the Fig. [I] for the input coherent field with 
the real amplitude ((f = 0). The estimated phase shift 
is always localized around the true value, but in general 
the phase estimation is biased. The accuracy may be 
easily assessed for an appropriately squeezed input state 
characterized by the wave function analogous to (Jill) 



P sq (x,0') = — exp{-e 

\/7T 



2rr 



^Hcosfl'] 2 }, (16) 



r being the squeezed parameter. The value r = corre- 
sponds to the coherent state. 

The phase information is the sharpest near the point 
9' = 7r/2. The resolution may be evaluated as 



Assuming further the optimum partition of squeezed 
state energy N = \a\ 2 + sinh 2 r as e r /2 ta \a\ w y/N/2, 
the optimum resolution corresponds to the well k nown 
ultimate limit of squeezed state interferometry l/V nN 2 . 
No optimization is necessary in the case of coherent 



state yielding the limit of coherent state interferometry 
1/y/nN. This statistical analysis well corresponds to the 
semiclassical (linear) approximation, when the phase res- 
olution is predicted by the intrinsic fluctuations of the 
signal 



A9' = AX(6') 



AX(9')) 
1 d9' 



Using the distribution (|16|) we conclude (X(9')) 
V2\a\ cos 9', AX (9') = e~7\/2 and therefore 

A6>CX : rp— — 77 . 



a sm( 



Semiclassical treatment represents good estimation in the 
regime of the best resolution, nevertheless it failures at 
the points close to 9' = 0. Moreover, here also the Fisher 
information tends to zero, since the quadratic term in the 
relative entropy (|TJ) disappears. The necessary assump- 
tions concerning the existence of the Fisher information 
are not fulfilled and, for example, the Cramer-Rao bound 
is not valid Ml. Nevertheless, the analysis of the phase 
distribution (O) shows that the phase resolution reads 



A0| / =o (X yj A4>\$>=n/2, 



yielding considerably worse phase resolution at this 
point. 

The block diagram of the phase detection based on the 
maximum likelihood estimation is sketched in the Fig. |^. 
Assuming the phases of the local oscillator and input co- 
herent states as 8 and cp, the balanced homodyne detec- 
tion measures the statistics of the quadrature operator 
(||) at the point 9 + ir/2. The phase difference may be es- 
timated using the likelihood function of measured phase 
sensitive data yielding the conditional phase distribution 



Pi 



homod 



oc exp[— 2n|a| 2 (sin0 — sin^')" 



(17) 



The predicted phase resolution is as in the Fig. |J, but 
shifted by the value n/2 in both the true and inferred 
phases. The best resolution is then achieved if 9 — ip = 0. 
For coherent input field the phase distribution at this 
point reads 



P 



0) oc exp(— 2n\a\ 2 sin 2 



Assuming further the total energy needed for such a real- 
ization of multiple measurement as n|a| 2 , the phase dis- 
tribution may be compared with the proposal of Vogel 
and Schleich Q. The relations ^ and ( pi] ) at the phase 
9 + n/2 yield in this particular case the phase distribution 
predicted by the maximum likelihood estimation. Never- 
theless, it need not be necessarily so in the general case of 
detection of an arbitrary phase shift, using an arbitrary 
input state or using an arbitrary phase sensitive variable. 
Estimation theory therefore naturally extends the "phase 
without phase concept" . 
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We demonstrated that any phase-sensitive measure- 
ment may serve for statistical prediction of phase shift. 
The content of phase information may be evaluated using 
the relative entropy of the phase in dependence on the 
observable probabilities only. The resolution predicted 
by the Fisher information is achieved if it exists. The 
proposed method based on the maximum likelihood esti- 
mation uses the information accumulated in the process 
of multiple measurement in the optimum way. This treat- 
ment better suits to the experimental conditions than 
sophisticated phase concepts. Particularly, the phase 
distribution obtained in the process of measurement is 
rather associated with nonlinear functionals (likelihood 
functional, relative entropy) than with the linear ones 
such as the distribution functions on the phase space 
are (Wigner function, Q-function). Moreover, since the 
phase distribution in realistic experiments is biased and 
phase shift dependent, the detailed statistical analysis 
free of any a priori restricting assumptions is always nec- 
essary. 
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FIG. 1. Phase distribution as function of inferred phase 
shift (f> in dependence on the true phase shift 6 for coherent 
input with total energy n\a\ 2 = 100. 

FIG. 2. Scheme of homodyne detection used for phase 
difference measurement. 
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